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Abstract: We formulate a prescription for computing Minkowski-space correlators 
from AdS/CFT correspondence. This prescription is shown to give the correct retarded 
propagators at zero temperature in four dimensions, as well as at finite temperature 
in the two-dimensional conformal field theory dual to the BTZ black hole. Using the 
prescription, we calculate the Chern-Simons diffusion constant of the finite-temperature 
M = 4 supersymmetric Yang-Mills theory in the strong coupling limit. We explain why 
the quasinormal frequencies of the asymptotically AdS background correspond to the 
poles of the retarded Green's function of the boundary conformal field theory. 
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1. Introduction 



The main prescription of the anti-de Sitter /conformal field theory (AdS/CFT) corre- 
spondence which allows one to compute the correlators in the boundary CFT, was 
originally formulated for Euclidean signature P, |J, and has been successfully used ever 
since. Working in Euclidean space is a common and convenient practice which usually 
does not lead to any restrictions since the results obtained in Euclidean AdS/CFT can 
then be analytically continued to Minkowski space, if desired. 

In many cases, however, the ability to extract the Lorentzian-signature AdS/CFT 
results directly from gravity is indispensable. Many interesting properties of gauge 
theories at finite temperature and density, most notably the response of the thermal 
ensemble to small perturbations that drive it out of equilibrium, can only be learned 
from real-time Green's functions. Even in the absence of a gravity dual of QCD, the 
AdS / CFT correspondence may provide helpful insights into properties of thermal gauge 
theories at strong coupling. 

In principle, one may try to avoid the Minkowski formulation of the AdS/CFT 
correspondence altogether by working only with the Euclidean version and using ana- 
lytic properties of the Green's functions to find the real-time propagators. In practice, 
however, such a method is not useful, because analytic continuation to the Minkowski 
space is possible only when the Euclidean correlators are exactly known for all Matsub- 
ara frequencies. Since gravity calculations in non-extremal backgrounds usually involve 
approximations of some sort (such as the high- and low-temperature expansions used 
in ref. [Q), exact expressions involving all Matsubara frequencies are normally be- 
yond reach. It is therefore crucial to have a prescription allowing one to compute the 
Minkowski correlators directly from gravity. 

Subtleties of the Lorentzian-signature AdS/CFT correspondence are well known 
P, EL U\. In this paper we shall not try to give a formulation of the Minkowski AdS/CFT 
in a form as general as the Euclidean version. Our goal is to formulate a working recipe, 
whose justification and full understanding would hopefully emerge in future research. 
The simple prescription that will be given is sufficient for the computation of the 
retarded propagator, and hence all other two-point correlators, but does not allow one 
to calculate higher-point correlation functions. We defer the treatment of higher-point 
correlators to future work. 

The paper is structured as follows. In section |2| we review the properties of the 
Euclidean and Minkowski thermal Green's functions. We formulate our recipe in sec- 
tion § and illustrate it on a simplest zero-temperature example. We compute the 
Minkowski-space correlators of the operator dual to the dilaton in AdS 5, and show that 
they agree with the analytic continuation of the classic Euclidean result. In section |] 
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another nontrivial check of the prescription is performed: we compute the retarded 
Green's functions in the two-dimensional CFT dual to the Banados-Teitelboim-Zanelli 
(BTZ) black hole background. We show that the result coincides with the analytic 
continuation from Euclidean space. We also find that the poles of the retarded Green's 
functions correspond precisely to the quasinormal frequencies of the general BTZ black 
hole, and we explain why this is to be expected in general in the framework of our 
recipe. Section [| is devoted the computation of the Chern-Simons diffusion rate in the 
the strongly coupled M = 4 supersymmetric Yang-Mills (SYM) theory. The result can- 
not be tested independently, and is considered as a prediction. Section || contains the 
conclusions and outlook. The Appendices are devoted to various calculations outside 
the main line of the paper. 

The prescription developed in this paper is further applied to the hydrodynamic 
regime of the M = 4 SYM theory in ref. JTO . 



2. Minkowski thermal correlators 

Since the main topics of this paper is the computation of thermal Green's functions 
from gravity, let us review the definition and general properties of different thermal 
Green's functions. Consider a quantum field theory at finite temperature. Let O be an 
arbitrary local operator, which we assume, for definiteness, to be bosonic (the fermionic 
case is completely analogous). In Minkowski space, the retarded propagator for O is 
defined by 

G R (k) = -i I d 4 x e~ ik - x 9{t) ([6(x), 6(0)}) (2.1) 



(we use the — h++ metric convention, and t = x°). The advanced propagator is defined 
similarly, 

G A (k)=i [d 4 xe- lkx 6(-t)([d(x),d(0)}} . (2.2) 



From the definitions can see that G R {k)* = G R (—k) = G A {k). If the system is P- 
invariant, then KeG RA are even functions of cu = k° and ImG R,A are odd functions of 
us. Instead of G R and G A , one can consider other correlation functions of O in thermal 
equilibrium. One example is the symmetrized Wightman function, 

G(k) = d 4 xe- ik - x (d(x)6(0) + 6(0)6{x)} . (2.3) 

All other correlation functions can be expressed via G R , G A and G. For example, the 
Feynman propagator is 

G F (k) =-ij d A x e~ ik - x (TO(x)6(0)) = ~ [G R (k) + G A (k)] - iG(k) . (2.4) 
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From the spectral representations of G R and G, one can relate G(k) with the 
imaginary part of the retarded propagator, 

G(k) = - coth ^ Im G R (k) . (2.5) 

Therefore, if G R (k) is known, all other correlators are easy to compute. In particular, 

G F (k) = Re G R {k) + i coth ^- Im G R {k) . (2.6) 



At zero temperature (T — > 0), eq. (|2.6| ) reduces to a simple formula, 



G F {k) = ReG R (k) + isgnu>lmG R (k) , T = 0. (2.7) 

Taking the limit uo — * in eq. ( |2.5| ), one obtains another useful formula, 

2T 8 
G(0, k) = - lim — Im G R (k) = 2iT—G R (to, k) . (2.8) 

The multitude of Minkowski-space Green's functions is in contrast with the situa- 
tion in Euclidean space, where one normally deals only with the Matsubara propagator, 

G E (k E ) = J d 4 x E e- lkE - XE (T E d(x E )d(0)) . (2.9) 

Here T E denotes Euclidean time ordering. The Matsubara propagator is defined only 
at discrete values of the frequency uj e . For bosonic O these Matsubara frequencies are 
multiples of 1i\T . 

The Euclidean and Minkowski propagators are closely related. The retarded prop- 
agator G R (k), as a function of u, can be analytically continued to the whole upper half 
plane and, moreover, at complex values of u equal to 2niTn, reduces to the Euclidean 
propagator, 

G R (2niTn, k) = -G E (2nTn, k) . (2.10) 

Analogously, the advanced propagator, analytically continued to the lower half plane, 
is equal to the Matsubara propagator at the points uj = —2iTiTn, 

G A {-2mTn, k) = -G E (-2ixTn, k) . (2.11) 

In particular, by putting n = 0, one finds that G R (0, k) = G A (0, k) = —G E (0, k). 
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3. Minkowski AdS/CFT prescription 



3.1 Difficulties with Minkowski AdS/CFT 

Let us first recall the formulation of the AdS/CFT correspondence in Euclidean space. 
For definiteness, we shall talk about the correspondence between the strongly coupled 
M = 4 SYM theory and classical (super)gravity on AdSs x S 5 . The Euclidean version 
of the metric of the latter has the form 

ds 2 = ^(dr 2 + dx 2 + dz 2 ) + R 2 dn 2 (3.1) 

z 1 

which is a solution to the Einstein equations. In the AdS/CFT correspondence, the 
four- dimensional quantum field theory lives on the boundary of the AdSs space at z = 0. 
Suppose that a bulk field (f> is coupled to an operator O on the boundary in such a way 
that the interaction Lagrangian is (f)0. In this case, the AdS/CFT correspondence is 
formally stated as the equality 

( e / 8M ^) = e -s«aM j (3.2) 

where the left-hand side is the generating functional for correlators of O in the boundary 
field theory (i.e., Af = 4 SU(N) SYM theory at large N and large 't Hooft coupling 
g 2 N), and the exponent on the right-hand side is the action of the classical solution to 
the equation of motion for in the bulk metric with the boundary condition 0| z= o = 4>o- 
The metric corresponds to the zero-temperature field theory. To compute the 
Matsubara correlator at finite temperature, one has to replace it by a non-extremal 
one, 

ds 2 = % { f(z)dr 2 + dx 2 + ^) + R 2 dQ 2 5 , (3.3) 
z 2 \ f[z)J 

where f(z) = 1 — z 4 /z H and zh = i^T)" 1 , and T is the Hawking temperature. The 
Euclidean time coordinate r is periodic, r ~ r + T _1 , and z runs between and zh- 

One can try to formally write the Minkowski version of the AdS / CFT as the equiv- 
alence 

(J'Ssm^Ok _ e iS cl [<p] _ 

However, an immediate problem arises. In the Euclidean version, the classical solution 
is uniquely determined by its value <p at the boundary z = and the requirement 
of regularity at the horizon z = Zh- Correspondingly, the Euclidean correlator ob- 
tained by using the correspondence is unique. In contrast, in the Minkowski space, 
the requirement of regularity at the horizon is insufficient; to select a solution one 
needs a more refined boundary condition there. This problem is well known ||, and 
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is thought to reflect the multitude of real-time Green's functions (Feynman, retarded, 
advanced) in finite-temperature field theory. Previous discussion of the Lorentzian 
signature AdS/CFT appears also in refs. |6|, [7], [|, 0. 

One boundary condition at the horizon stands out from the physical point of view. 
It is the incoming-wave boundary condition, where waves can only travel to the region 
inside the horizon of the black branes, but cannot be emitted from there. We may 
suspect that this boundary condition corresponds to the retarded Green's function, 
while the outgoing-wave boundary condition gives rise to the advanced Green's function. 
However, even after fixing the boundary condition at the horizon, it is still problematic 
to get eq. ( |3.4|) to work. 

To see where the problem lies, first we notice that the AdS part of the metric ( p.3|) 
is of the form 



Let us consider fluctuations of a scalar (e.g., the dilaton or the axion) on this back- 
ground. The action of the scalar reads 



where K is a normalization constant (for the dilaton K = —tc 3 R 5 /4k>1 , k,iq is the 
ten-dimensional gravitational constant) and m is the scalar mass. The limit of the 
z-integration is between the boundary zb and the horizon zh (for the metric Q3.3|), 



ds 2 = g zz dz 2 + g^ LV {z)dx^dx 1 ' . 



(3.5) 




(3.6) 



Z B = 0). 



The linearized field equation for 0, 




(3.7) 



has to be solved with fixed boundary condition at z B . The solution is 




(3.8) 



where <po(k) is determined by the boundary condition, 




(3.9) 



and fk(z) is the solution to the mode equation, 




(3.10) 



6 



with unit boundary value at the boundary, f k {zB) = 1 and satisfying the incoming- wave 
boundary condition at z = zjj. 

The action on shell (i.e., for classical solutions) reduces to the surface terms 



d A k 



S = 77TTiM-k)Hk,z)Mk) , (3-11 



(27T) 



Z = Z H 



Z=Z B 



where 

F(k, z) = K^g zz f_ k (z)d z f k (z) . (3.12) 



Now, if one try to make the identification (|3.4j) , the Green function can be obtained 
by taking second derivative of the classical action with respect to the boundary value 
of (ft. From eq. ( 3.11|) , this Green's function is 



ZII 



G(k) = -F(k,z) -T{-k,z) . (3.13) 



ZD 



ZII 



ZD 



One can see, however, that this quantity is completely real, and hence cannot be a 
candidate to the retarded Green's function, which is in general complex. Indeed, by 
noticing that f k (z) = f- k (z) (this is because f k (z) is also a solution to eq. ( |3.10|) and 



also satisfies the incoming- wave boundary condition), one sees that the imaginary part 
of T{k, z) is proportional a conserved flux, 

ImF(k, z) = %zV=gg" [r k d z f k - f k d z f* k ] , (3.14) 

and thus \m.J r {k,z) is independent of the radial coordinate, dzlmj 7 = 0. Therefore, 
in each term in the right-hand side of eq. ( |3.13| ), the imaginary part at the horizon 
z = zh completely cancels out the imaginary part at the boundary z = zb- 

We should clarify one issue with the boundary term at the horizon. It is commonly 
thought that this term is oscillating and hence averaged out to zero. However, this is 
not the case if the incoming-wave boundary condition is consistently maintained. We 
shall see that explicitly below in section |3.3| . 

We can try to avoid this problem by throwing away the contribution from the 
horizon, keeping only the boundary term at z = Zb- However, now the imaginary parts 
of the two terms in eq. ( |3.13|) cancel each other: from the reality of the field equation 
one can show that J-(—k,z) = T*{k,z\ Therefore, the resulting G(k) is still real. It 
seems that by differentiating the classical action one cannot get the retarded Green's 
function, which is, in general, complex. 
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3.2 The Minkowski prescription 

We circumvent the difficulties mentioned above by putting forward the following con- 
jecture 

G R (k) = —2T{k, z) . (3.15) 

ZB 

which seems rather natural but, for reasons explained above, does not follow strictly 
from an identity of the type (|3.4j) . The justification for eq. ([3.15|) is that it works in all 



cases where independent verification is possible. 

Our prescription for computing the retarded (advanced) Green's functions is thus 
formulated as follows: 

1. Find a solution to the mode equation ( |3.10| ) with the following properties: 

• The solution equals to 1 at the boundary z = z B \ 

• For timelike momenta, the solution has an asymptotic expression corresponding 
to the incoming (outgoing) wave at the horizon. For spacelike momenta, the solution 
is regular at the horizon. 

2. The retarded Green's function is then given by G = — 2J-"qm, where T is given 
in eq. ( |3.12| ). Only the contribution from the boundary has to be taken. Surface terms 
coming from the horizon or, more generally, from the "infrared" part of the background 
geometry (corresponding to the "position of the branes") must be dropped. This part 
of the metric influences the correlators only through the boundary condition imposed 
on the bulk field <j). 

Since the imaginary part of J-(k, z) is independent of the radial coordinate, Im G R,A 
can be computed by evaluating lmJ-(k, z) at any convenient value of z. In particular, 
it can be computed at the horizon. 

In order to verify that the prescription indeed gives the retarded propagator G R , 
one should compute G R in theories where it is known from other methods. Below we 
perform a check for zero-temperature field theory; further checks are done in section (|J) 



and ref. 10 



3.3 Example: Zero-temperature field theory 

As an example, we use the prescription ( |3.15| ) to compute the retarded (advanced) 
two-point function of the composite operators O = \F 2 at zero temperature. In this 
case the action (|3.6|) is that for a minimally coupled massless scalar on the background 



(|3.1|) . The horizon is at zu = oo and the boundary at zb — e — ► 0. The mode equation 



reads 

f',!(z)-- z f k (z)-k 2 f k (z) = 0. (3.16) 
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For spacelike momenta, k 2 > 0, the calculation is identical to the one for the Euclidean 
case (see also appendix [A]), the only difference being the extra minus sign in front 
of the Lorentzian signature action. We obtain 

N 2 k 4 



G R (k) 



64vT 2 



Ink 2 



k 2 >0. 



(3.17) 



For timelike momenta, we introduce q = \J—k 2 . The solution to eq. ( |3.16| ) satisfying 
boundary conditions outlined above is given by 

' z 2 H^\qz) 



fk(z) 



e 2 H?\qe) 
z 2 Hf\qz) 



u > 0, 



io < 0. 



(3.18) 



Notice that /_ 



V^Hf\qe) 

/£. Computing (|3.12| ) and using eq. ( |3.15|) , we get 

N 2 k A 



G R {k) 



647T' 



-(In A; 



nx sgn LU ) 



Combining eqs. ( |3.17 ) and ( 3.19 ), we obtain the retarded Green's function, 

N 2 k A 



G (k) 



64tt 2 



(In \k | — m6(—k ) sgncu) . 



(3.19) 



(3.20) 



Using the asymptotics of the Hankel functions for z — > oo, one can check that !F(k, z) is 
not zero at large z — > oo. Moreover, it is purely imaginary in this limit, and is equal to 
iiV 2 fc 4 sgnu;/1287r = ImjF(A;,e), as it should be in view of the flux conservation fl3.14|) . 
Thus, the imaginary part of the Green's function can be computed independently using 
the asymptotics of the solution at the horizon. 

By using the relation between the Feynman and retarded propagators at zero tem- 
perature, eq. Q2.71) , we find 

N 2 k 4 



G*{k) 



64vr 2 



(ln|fc 2 | -in9(-k 2 )) 



which can be obtained from the Euclidean correlator, 



G E {h 



(3.21) 



(3.22) 



by a Wick rotation. Thus, our prescription reproduces the correct answer for the 
retarded Green's function at zero temperature. 

The Euclidean and Minkowski correlators corresponding to massive scalar modes 
in AdS 5 are computed in appendix 
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4. Retarded correlators in 2d CFT and the BTZ black hole 



In this section, we use the prescription ( p.15 ) in the general BTZ black hole back- 
ground to compute the retarded Green's functions of the dual two-dimensional CFT. 
The result is then analytically continued to complex frequencies, and is shown to repro- 
duce the Matsubara correlators known from field theory. The calculation is technically 
quite involved, but is nevertheless instructive since it is the first nontrivial check of our 
prescription at finite temperature. As a by-product, we demonstrate that the gravi- 
tational quasinormal frequencies in the asymptotically AdS background correspond to 
the poles of the retarded correlator in the dual CFT on the boundary, and thus provide 
a quantitative description of the system's return to thermal equilibrium. 



4.1 Gravity calculations 



The non-extremal BTZ black hole |TT], [L2| is a solution of the 2 + 1-dimensional Einstein 
equations with a cosmo logical constant A = — l/l 2 . The metric is given by 



p P 2 



i P 2 



i 2 P 2 



+ 



(p 2 - pD(p 2 - p 2 -) 



dp 2 , (4.1) 



where <p is periodic with the period 2tc and p± correspond to the positions of the 
inner /outer horizons. Introducing new variables p, x±, with 



x 



p 2 + cosh p — p 2 _ sinh p 
± P±t 



(4.2) 



we can write the metric as 

ds 2 = l 2 dp 2 — sinh 2 p(dx + ) 2 + cosh 2 p(dx~) 2 



(4.3) 



In terms of z — tanh 2 p, the action of a massive scalar in the BTZ background reads 



S = —7] / dz dx + dx 



I 



2 1- z) 



4z(l - Z) 2 2 l-Z 2 

-{d z <t>) (9+0) 



f- 



+ (1 - z) (d_0) 2 + m 2 <j> 2 



(4.4) 



where r] is the normalization constant. 

The mode functions fk + k-{ z ) are defined so that 



(j)(x d 



dk + dk- 

"W" 



-i(kj,.x + +k-x ) 



fk+,k-(z)M k +, k -) 



(4.5) 
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and satisfy the equation 

d 2 zz fk+,k- + -d z fk + ,k„ + 



2 U2 



l 2 k 



2l„2 



Vk 



m 2 l 



2/2 



Az 2 {l-z) 4z(l-z) Az(l-z) 2 



fc_i_,fc_ 



(4.6) 



On shell, the action ( O ) reduces to the boundary terms 
'dkj.dk- , 



S 



k+, -kJ) [F{k + , z B ) - F{k +1 z H )\ 4>o(k +} kJ) , (4.7) 



where zb = 1 corresponds to the boundary at infinity (p = oo), = is the location 
of the outer horizon (p = p + ). The function T is given by 



2z 



T{k+,k_,z) = -V^fk+,k. d zfk+,k 



(4. 



According to the Minkowski AdS/CFT prescription, the retarded/ advanced Green's 
functions are determined by 



G R ' A (k + , k- 



-2[lim F{k + , k-,z)] 



(4.9) 



where the solution fk + ,k-(z) is normalized to 1 at z — 1 — e and represents a purely 
incoming (outgoing) wave at the horizon z — 0. Square brackets indicate that the 
contact terms are ignored. The imaginary part of G R,A , due to the flux conservation, 
can be computed either at the boundary or at the horizon, 



\mG RA (k+,kJ) = -2[lim lmj r (k+, fc_, z)} = - 

Z^Z B 

The solution to eq. ( |4.6|) can be written as 

fk + ^{z) = z a {l-zfe-\l-e)- a 

where 



•2[lim lmF(k + ,k_,z)} . (4.10) 

z-^zu 



aFi (o,6; c; z) 
3F1 (a, b; c; 1 - e) 



a, 6 



/(fc+ =f fc_; 

2i 



fi. 



2a 



(4.11) 
(4.12) 

(4.13) 



and the indices a, (3 are given by 

a± = ±ilk+/2, /3± = ^ (l ± VI + / 2 m 2 ) = ^ , 

a_ corresponds to the incoming wave at the horizon, a + to the outgoing one, and A± is 
the conformal weight of the boundary CFT operator. The range of A + is [1, 00), while 
the range of A_ is (0, 1], where zero corresponds to the unitarity bound in d = 3. Here 
we shall consider only the A + branch. Accordingly, 1 we choose (3 = (3- in eq. ( 4.11] ). 

1 In AdS/CFT, if a supergravity scalar behaves near the boundary as / ~ Ae d ~ A+ + Be A +, A is 
interpreted as the source of a dimension-A + operator. 
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4.1.1 The full retarded Green's function 



Near the boundary, i.e., at z — 1 — e, the derivative of the solution (|4.11|) takes the 
form 

,/ / Ni P , a62Fi(a + l,6 + l;c+l;l-e) 

h + ,kA z )\z=i-e — 1 ^tt — r — i \ > l 4 - i4 J 

+ ' e c^i (a, b; c; 1 - ej 

The retarded Green's function is then computed using the prescription ( |4.9| ). 



G R (k+,k. 



T 



afc aFi(a + l,6+l;c+l;l-e) 
— hm — - 

c e— o 2 F 1 (a, o; c; 1 - e) 



(4.15) 



Some care should be taken when computing the limit in ([4.15 ) for 1 — 2/3+ = c — a — b = 
—n, with n = 0, 1,2, i.e., for A + = 2(3+ ^ n + 1. In this case the hypergeometric 
function is degenerate and contains logarithmic terms. Accordingly, integer and non- 
integer A + need to be considered separately. 



4.1.2 Non-integer A + 

For A + = 2j3 + ^ n + 1, where n = 0, 1, 2, direct calculation gives 



Gr(p+,P- 



7r/r 2 (2/3+ - l)sin2?rA 
P+ P- 



r ha 



cosh , 

2T L 2T R 
+ % sin 2tt(3 + sinh 



Ph 



cos 27T/9+ cosh 
P- 



*P+ 
2ttT l 

P+ 



r /a 



2^77 



+ 



P- 



277, 2T, 



(4.16) 



277 2Ti 



where T ivR = (p + =)= p_)/27r, and the momenta p± are related to k± and tu, A; as follows 

u + k 



p + = nT L (k + + k-)l= , 



ttT r (k+ - k-) I 



(4.17) 



4.1.3 Integer A + 

In the case of integer conformal dimension, one has to further distinguish between A + 
being an odd or even integer. In the latter case, i.e., for hi = h,R = f3+ = A + /2 = n+1, 



n = 0, 1, 2, we get 



G R (p+,P- 



2d 



l7r 2 T 2 (2n+ 1) 
sinh — ^- sinh ■ ^ 



r 1 +n + 



*P+ 
27rTr 



r 1 + n + 



2nT, 



R 



2T T 



2T< 



R 



ip ( 1 + n — 



2ltTr 



+ ip [1 + n- 



2ttT, 



R 



(4.18) 
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In particular, the following simple result is obtained for Hl — h R — 1 (A + = 2): 



VP+P- 
' n 2 lT L T R 

The imaginary part of ( |4.19|) , 



V I- 



2nT T 



+ V 1- 



2ixT, 



R 



(4.19) 



VP+P- 



2ixlT L T R 



coth + coth ^ 



277. 



277 



(4.20) 



coincides with the appropriate limit of eq. ((4.24Q below. For odd integer values of A + , 
i.e., for h L = h R = (3+ = A + /2 = n + 1/2, n — 1, 2, we get 



Gr(p+,p- 



4r]e 



2n-l 



ln 2 T 2 (2n) 



T[ 1 - + n + ^±- 
1 2 2nT L 



cosh cosh ■ ^ 



2T L 2T R 
Finally, for (3+ = 1/2 (A+ = 1) we have 



F(i + n+^- 
1 2 2ttT r 



— — + w \ — Vn — 

2nT L ) v \2 2nT R 



(4.21) 



G R (p + ,p- 



4r] 
le\n 2 e 



' 2 2ttT l 



1 2 2ixT R 



(4.22) 



The imaginary part of the obtained Green's functions coincides with eq. ( (4.24Q below. 
More crucially, one can check that the only singularities of ( [4.16| ), ( |4.18| ), ( |4.21| ) and 
( |4.22| ) are simple poles located in the lower half-plane of complex u at 



CO 



W = k-vlirT L (h L + n), n = 0,1,2,.. 
-k — i4irT R (h R + n) , n — 0, 1, 2, 



(A) 



(4.23a) 
(4.23b) 



No singularity is located in the upper half plane, as expected for the retarded Green's 
function. 



4.1.4 A note about the imaginary part 

Notice that if we are interested only in the imaginary part of the retarded propagator, 
the calculations are much simpler because it can be done by using the solution (|4.11|) 
with a = a_, (3 = j3_ and taking the limit z — > Zh in eq. ( (4.101 ). The result reads 



Im G R = - C (e, h LR ) sinh ( + V ~ 



277 2T< 



R 



r h L + 



2nT r 



r h R + 



%P- 



2nT, 



R 



(4.24) 
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where = Hr = @+. The normalization constant is 

C(e,h LR ) = <*ir(2ft L -l)r(2/*-l)' A+>1 ' (4 . 25) 
2r ? a _i 

, i 2 ' ZA + - 1. 

7r / e m e 



This result has been known before from the absorption calculations [13 



4.2 Comparison with the two-dimensional finite-temperature CFT 

Now let us show that the retarded Green's function found from gravity coincides with 
the field-theory result. One starts from the coordinate-space expression for the Eu- 
clidean finite-temperature correlation function of the local operator O with conformal 



dimensions (h, h) 



G V, .) = ( („, ,) O (0, 0)) = ^ ) , (4,6) 



where Co is a normalization constant and w = x + it, with x running from — oo to 
+oo, and r being periodic with the period 1/T. 

To compare with the results obtained from gravity, we need to compute the Fourier 
transform of (|4.26|) for h — h, obtain the momentum-space Matsubara correlator, 



and compare the result with the retarded Green's functions obtained from gravity in 
section |4.1.1| . We make the comparison only for Tl = Tr and for the integer values of 
A = 2h = 2h. We need to compute the following integral, 

l/T 

G E {u E = 2vmT, k)= fdr [ dx e iuJET e ikx G e (t, x) . (4.27) 



It is convenient to write G (r,x) in the form 

G£(T ' X) = (cosh 2nTr- cos 2ttTt)* ' (4 ' 28) 

We shall first take the integral over r and then over x. Introducing z = exp (— i27rTr) 
and r = \x\, the integral becomes 

oo 

;-l) 2A 2(27rT) 2A C £ dz ^gA—n—i 

2+ 



G (ue, fc) = r drcctkr f - j—-—^ , (4.29) 



|*|=l 
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where z± = exp (±27rTr). The integral over z is computed by evaluating the residues. 
The result reads 



G E {u E ,k) 



2(2irT) 2A C T(A-n) 
f r(A) 



A-l 

£ 

m=0 



r(2A — m — V 



m!T(A — m)T(A — m — n 



■J n , m (A;), (4.30) 



where J n>m {k) are the following integrals that can be computed explicitly |I5 



Jn,m\k) 



dr 



e 2nTr(A+n) CQS fcf 



\m-2A 



r( s ) 



\2A-m-l 



8ttT 



r(z) 



(4.31) 



and s = 2 — 2A + m, and 



3A n zA; 

2 hmH — 

2 2 4ttT 



(4.32) 



In general J n , m (k) diverges, and the result on the right-hand side of eq. (|4.31]) should 
be understood as a limit of the regularized expression; after taking the limit, the (di- 
vergent) contact terms should be ignored. 

As an example, consider the simplest case A = 1. Taking the limit s — ► in 
eq. fl4.31p , we get 

1 



Jnfi{k) 



8ttT 



Therefore, we obtain 

G E {2nnT,k) = -7i C 



■0 



l + n 



ik 



1 + n ik 
— ;— + 



AttT 



(4.33) 



(4.34) 



Identifying the prefactors as Arj/leln 2 e = ttCq one observes that the retarded Green's 
function ( f4.22|) taken at uo = 2iriTn equals (with a minus sign) to the expression given 
by eq. ( |4.34| ), as expected on general grounds in eq. ( |2.10| ), and therefore the field theory 
and gravity results are consistent. 

For A = 2, eqs. QOD|) and ( CT ) give 



G E (2%nT, k) = nC { 



a- 



(2vrTn) 2 + k 2 



. n ik 

ib H 

1 2 AtxT 



. n ik 

+ ib H 1 

1 2 4ttT 



(4.35) 

Upon identification tc'Cq = v/ITrT l , ( |4.35|) coincides with — G r {uj = 2rrTn, k), where 
G R was computed in eq. (|4.19|) . One can repeat this procedure for any integer A. So 
we have checked that our prescription does give the correct retarded Green's functions 
in two-dimensional CFT. 
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4.3 Quasinormal modes and singularities of G R 

The quasinormal modes for black holes in asymptotically AdS spacetimes are defined as 
the solutions to the wave equation obeying the "incoming wave" boundary condition at 
the horizon and the vanishing Dirichlet condition at the boundary. 2 It was suggested 



in ref. |18| that the associated quasinormal frequencies are related to the process of 
thermalization in the dual strongly coupled CFT. 3 The inverse relaxation time is given 
by the imaginary part of the lowest quasinormal frequency, which has been computed 



numerically in ref. |18| for AdS black holes in various dimensions. 

Due to the simplicity of the BTZ background, the quasinormal modes in this case 
can be determined analytically p5[ |27j. It was first noticed in ref. p8| that the quasi- 



normal frequencies found in ref. |25| correspond to the poles of the retarded Green's 
function describing the relaxation of the boundary CFT to thermal equilibrium, as in 
the linear response theory. 

Our calculation of the retarded Green's function confirms this conjecture explicitly: 
the poles of G R (eqs. ( |4.23| )) are precisely the quasinormal frequencies of the BTZ black 
hole |§ . 



In fact, given the prescription (|3.15| ), it is not difficult to see that this should be true 
in general. Indeed, suppose <f>k(z) = A(k)4>\{z) + B{k)(f)2(z) is the (super)gravity field 
whose boundary value 0o couples to an operator O in the dual theory in Minkowski 
space, and suppose further that obeys the incoming wave boundary condition at 
the horizon. The normalized solution used to compute the retarded correlator has 
the form ipu(z) = o (^)0fc(^)/0fc( e )- Near the boundary z = e — > one has ipk( z ) = 
A(k) z A - +■■■ + B(k) z A + + ■■■. Then 

G R (k) ~ VV = (A+ - A_)^|e A +- A — 1 + contact terms + 0(e) . (4.36) 

Zeros of A(k) correspond to singularities of the retarded correlator. But A(k) = is 
precisely the vanishing Dirichlet boundary condition at z = which defines the quasi- 
normal modes in asymptotically AdS spacetimes (recall that the boundary condition 
at the horizon is the same for the quasinormal modes and for the solution from which 
G R is determined). Thus, quite generically, quasinormal frequencies are equal to the 
the resonance frequencies of the retarded Green's function of the operator O. In the 
linear response theory, they characterize the time scale of the approach to thermal 
equilibrium in the boundary field theory perturbed by O. 



2 For reviews and references on quasinormal modes in asymptotically flat spacetimes sec Jig , ^7| . 
3 Similar ideas appear in ref. [L9[. Other works on quasinormal modes in asymptotically AdS 



spacetimes include refs. 
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5. Thermal M = 4 SYM theory and the Chern-Simons diffusion 
rate 



We now consider the four- dimensional Af = 4 SYM theory at finite temperature. The 



thermodynamics of this theory has been studied in detail by using its gravity dual |29|, 
p0|1 . Here, we compute a simple thermal correlation function by using the Minkowski 
AdS/CFT prescription. 

The quantity we shall be interested in is the Chern-Simons diffusion rate, 

r= J d'x (d(x)d(o)) , (5.i) 

where 

6 = \f; u f^ . (5.2) 

This quantity is interesting because it has been intensively studied at weak coupling. 
In the standard model, V determines the rate of anomalous baryon number violation 



at high temperatures. Its knowledge is important for electroweak baryogenesis pTJ . At 



weak coupling, the parametric behavior of V is B3|, |35 



T = const • (g YM N) 5 In - ir — T 4 . (5.3) 

9ym n 

The overall numerical constant is of nonperturbative nature but can be found by a 
lattice simulation of a classical field theory |36|, [37], [38| . No method has been developed 
to compute T at strong coupling. 

In the original AdS / CFT correspondence, O is coupled to the bulk axion, therefore 
the correlators of O can be found from gravity. Indeed, if one finds the retarded 
propagator G R for O, then, according to eq. (|2.8|), 

9ym\ i;„ t „ nR 



, . lim — ImG K (cj,0). (5.4) 

Therefore, the task of computing T is reduced to the calculation of the retarded Green's 
function of O at small momenta. 

The coupling of the axion (the Ramond-Ramond scalar) Co to O is determined 



from the three-brane Born-Infeld action [^, . The action for Co is the same as for a 



minimally coupled scalar, eq. (|3.6|). It is convenient to use the coordinate u defined as 



z 2 



u = -j- • (5.5) 



^2 
Z H 
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The mode equation has the form 



K- 



l + u 2 



:f k + 



or 



u(l-u 2 ) Jk ' u{l~u 2 ) 2 
where the following notations were introduced 



fk 



u(l — U 2 ) 



f k = o, 



\k\ 



vo 



(5.6) 



(5.7) 



2vrT ' 2nT 

Equation (|5.6| ) is the Heun differential equation |4l|] (second-order differential equation 
with four regular singular points, see |42], |43|). Global solutions to the Heun equation 
are unknown. However, for computing the Chern-Simons diffusion rate one needs to 
consider only the regime of low frequency and momentum, u, \k\ <C T. In this regime, 
the solution representing the incoming wave at the horizon is given by 



fc (n) = (l- M )-^ 2 F fe (n), 
where Fk(u) is regular at u — 1 and can be written as a series 

1 + u 



Fkiu) 



T +C 2 )ln. 



0(w 2 , W, q 4 ) • 



(5.8) 



(5.9) 



Computing the boundary term ( |3.12| ) near u = we obtain 



vr 2 iV 2 T 4 



^(k, u) 



ixo 



u 



+ 0(t» 2 , n>q 2 , q 4 



When k = 0, our prescription implies 

G R/A (oo,0) = T 
From eq. ( |5.4|) one then finds 



inN 2 uT 3 



1 + 



rp4 



(5.10) 



(5.11) 



(5.12) 



256vr 3 

which is our final result for the Chern-Simons diffusion rate at large 't Hooft coupling 
gy M N 3> 1. Note that (|5.12|) has the same large N behavior as the weak-coupling 
expression ( |5.3|) in the 't Hooft limit (that is, of order N°). The natural conjecture to 
make is that the Chern-Simon diffusion rate is proportional to T 4 with a proportionality 
coefficient dependent on the 't Hooft coupling, 



where fr{x) ~ x 5 In - for iCl and fr{x) = x 2 / (2567T 3 ) for x ^> 1. 



(5.13) 
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Notice that if we considered, instead of O = \FF the operators \F 2 , the answer 
would be the same. The same result is also valid for spatial components of the stress- 
energy tensor if k — 0, since the corresponding gravitational perturbation obeys the 
same equation ||40|| . For the latter case, we reproduce, from Kubo's formula, the value 
for the shear viscosity rj 



f A 2 T 3 that has been found in ref. jH 



The imaginary part of the retarded propagator G(u, k) can also be computed for 
very high frequencies and momenta [w, \ k\ 3> T). This is done in appendix [B]. 



6. Conclusions 

We have formulated a prescription that allows one to compute the Minkowski retarded 
Green's function from gravity. We have demonstrated that the prescription works 
for zero-temperature gauge theory, and for finite-temperature two-dimensional CFT. 
By using the prescription, we found the Chern-Simons diffusion rate for the M = 4 
SYM theory at finite-temperature. The prescription allows one to compute many other 
quantities which cannot be computed from Euclidean space. Two examples, the shear 
viscosity and the R-charge diffusion rate in thermal TV = 4 SYM theory, are presented 
in ref. |IU| . 

The prescription given here does not allow one to find higher-point correlation 
function, and hence it cannot be the final word on Minkowski gauge theory /gravity 
duality. Hopefully, it can be incorporated in some future more general scheme. 
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A. Zero-temperature correlation functions in 
Euclidean and Minkowski AdS/CFT 

A.l Euclidean signature 

To compute the Euclidean two-point function of a CFT operator O one uses the 
AdS/CFT equivalence 

( e /sM^°) = e" SEM , (A.l) 

where Se[4>} is the classical (super)gravity action on M and 4>o is the boundary value 
of the bulk field 0. At zero temperature M = AdSs x S 5 . These calculations were first 
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done in ref. and then appeared in the literature many times. We repeat them here 
to have an explicit comparison with the Lorentzian version. The (Euclidean) AdSs part 
of the metric is given by 

R 2 



ds d+1 = —^{dz + dx 



(A.2) 



where x are the coordinates on R A . The action of the massive dilaton fluctuation <fr is 

vr 3 /? 8 



S r 



Ak 2 



dzd xz 



(A.3) 



Using the Fourier representation, 

<j)(z,x) = 

we obtain 

oo 

tt 3 /? 8 f , [d 4 kd A k' S(k + k r 



d A k 



e lkx f k {z)Uk) , 



(A.4) 



S 



dz 



(27T) 



m 2 R 2 

d z fkd z fk' - kk'f k f k , H 5— /fc/jb' 



The equation of motion for fk{z) reads 



ft- -fx- 1^ 



m 2 i? 2 



fk = 



It has a general solution 



(z) = Az 2 I v (kz) + Bz 2 I^ v {kz) 



(A.5) 
(A.6) 

(A.7) 



where v = \/4 + m 2 R 2 (or equivalently one can consider I v and as two independent 
solutions). We also use the standard notation A = v + 2 for the conformal weight of 
the operator O. The solution regular at z = 00 and equal 1 at z = e is given by 



fk(z) 



z 2 K v {kz) 



e 2 K u (ke) ' 

On shell, the action reduces to the boundary terms, 



(A.8) 



4-K 2 



7T 3 i? 8 fd A kd A k' 



(27T) 



(2tt) V(fc + k')<f> (k)Mk 



r J'k'{z)d z f k {z) 



d 4 kd 4 k' 

(2vr) 8 



(A.9) 
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We note that J-{z, k, k') is a real function, and JF(oo, k, k') = 0. The two-point function 
is given by 



(6{k)6{k')) = z 



-2F(z, k, k') 



</>o=0 



(27r) 4 5 4 (fc + k') 



,yR s f k >(z)d z f k (z) 



(A.10) 



Direct calculation from eq. (|A.8|) gives 

n 3 R 8 



(0(k)0(k')) 



2k 



, e 2 ( A - d )(27r) 4 5 4 (A; + k')k 2u 2 1 ~ 2u ^^ / u) + 



10 



(A.ll) 



where dots denote terms analytic in and/or those vanishing in the e —>■ limit. 
Substituting re 10 = 2rr 5 / 2 R 4 /N || we get 



For integer A, the propagator is given by 

(_1)A N 2 



(A.12) 



(0(k)0(k>)) 



pA-4 

e 2(A-4) (27r) 4 5 4 (fc + ^^ ,,. ,.2 



hr, A = 2,3,4, 



(A-3)! 2 8vr 2 v y v y 2 2A " 5 

(A.13) 

where one sets (—1)! = 1. Specializing further to the massless case (A = 4), we have 



7V 2 

(0(k)0(k')) = -^-(2K) A 8\k + k'Wlnk 2 . 
A. 2 Lorentzian signature 

In the AdSs background the action describing the dilaton fluctuation reads 



(A. 14) 



S 



TT 3 R 8 



d x I dz 



dz 



(A.15) 



where 77^ = diag(— 1, 1, 1, 1) and R is the AdS§ radius. Writing the classical solution 

as 



[z, t, x) 



d 4 k 

the action reduces to the boundary term 

d 4 k 



-iiot+ik-x 



f k (z)(f) (k) 



S 



(2tt)< 



T{k, oo) - Tik, e) (f) (-k)(f) (k) 



(A.16) 



(A.17) 
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where 



F(k,z) 



ir 3 R 8 f- k d z f k (z) 
4«i z 3 



(A.18) 



The mode function fk{z) obeys the same equation (|A.6|) as in the Euclidean case. 
However, now the form of the solution depends on whether k is spacelike or timelike. 

A. 2.1 Spacelike momenta, k 2 > 

For spacelike momenta, the solution is real and the calculation is identical to that of 
the Euclidean case, except for the additional minus sign coming from the difference in 
definitions of Euclidean and Minkowski actions. The result is thus given by 



G R (u, k) 



N 2 (k 2 ) A - 2 T(3-A)e 2 ^ 



8vr 2 2 2A - 5 r(A - 2) 



(A.19) 



A. 2. 2 Timelike momenta, k 2 < 



Introducing q = \/ — k 2 = \k\, the solution to ([A. 6 ) satisfying the "incoming wave' : 
boundary condition at z = oo and normalized to 1 at z = e is 



4>k(z) 



z 2 Hj l \qz) 
e 2 Hi l \qe) 
z 2 H<?\qz) 
U 2 i#(ge) 

Here v = \J A + m 2 R 2 . Using the expansion 



io > 



u < 0. 



(A.20) 



z 2 H^ 2 H 



Sin Z/7T 



2 V Z 



n-i/ 



(l + 0(z 2 ))+e^ 



yV + 2 



2T(i/ + 1) 



(l + 0(, 2 )) 



(A.21) 

and the prescription (|3.15|) , one obtains the retarded Green's function for timelike 
momenta 



G R (k) 



N 2 {q 2 ) A - 2 T{3 - A)6 2 ( A ' 4 ) 
8vr 2 2 2A - 5 r(A - 2) 



[cos 7rA — i sin (ttA) sgn uj] . 



(A.22) 



Here A = 2 + z/ 7^ 2, 3, .... For integer values of A, by taking the appropriate limit in 
eq. ( [A.22Q , one has 



G R {k) 



iV 2 e 2(A-4) 



8vr 2 (A-3)! 2 2 2A - 5 
where one sets (—1)! = 1. 



(9 



2\A-2 



[in g 2 — Z7T sgn u\ , A 



2,3,4, 



(A.23) 
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Combining the results for the time- and spacelike momenta, one can write 

k 2 > 0, 
cos7rA , k 2 < 0, 



ReG R {k) 



iV 2 |A; 2 | A - 2 r(3- A)e 2 ( A " 4 ) 

87r 2 2 2A - 5 r(A - 2) 
iV 2 |A; 2 | A - 2 r(3- A)e 2 ( A - 4 ) 



(A.24a) 



8vr 2 2 2A - 5 r(A - 2) 



7\r2|jt.2|A-2TYo _ A\,2(A-4) 

Im G R (k) = 8 J 22 2A- 5r(A _ 2) sin (* A ) e (~ k2 ) s S n " • ( A - 24b ) 

For the integer values of A, we have, correspondingly, 

N 2 2(A-4) 

ReG (fc) = - — oXI9o9A _J fc 2 | A - 2 In |fc 2 | , A = 2, 3, 4, . . . (A.25a) 



8tt 2 (A-3)! 2 2 2A - 



jyr2 2(A-4) 

ImG (&) = ~ 87r(A _ 3) ,2 2 2A-5 l fc2 | *(-* 2 ) s S n w > A = 2, 3, 4, . . . (A.25b) 
In the massless case (A = 4), we get 

N 2 k A 

G R (k) = T (ln\k 2 \ -m9(-k 2 )sgiiu). (A.26) 

647T 2 

B. Retarded Green's function at high momenta 

In this appendix we consider the retarded Green's function of O = \F 2 (the same 
formulas apply for O = \FF) in the regime u, \k\ ^> T. Since the temperature now is 
much smaller than the momentum scales, one should expect the result to be close to the 
zero-temperature result (|3.20|) , with the temperature providing only a small correction. 



This can be shown explicitly by using the Langer-Olver method 0, 45]. However, one 
qualitative effect occurs at finite temperature: G R (k) acquires a small imaginary part 
for spacelike fc's (recall that at zero temperature, G R (k) is real if k 2 > 0). This happens 
for the same reason as the one causing the Landau damping in plasma: excitations in 
the plasma can absorb a spacelike momentum. We shall estimate lmG R (k) for large 
spacelike fc's. 

By a change of variables from to W(u) = u~ l l 2 \/u 2 — 1 fk{u), the mode 

equation (|5.6|) can be written in the form of the Schrodinger equation, 



W"(u) = (vo 2 F(u) + G(u)) W(u) , (B.l) 

where 

\-s\l-u 2 ) n 4 + 6n 2 -3 
F W = T\ > G(u) = —-r~2h 2^2' B ' 2 
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and s = \k\/tj. In the regime, vo = (uR) 2 /2p = uj/2tiT ^> 1, s ~ 1, the term 
proportional to ix> 2 dominates in the potential, and the solution can be found by the 
WKB approximation. 

For s > 1 (corresponding to timelike fc's), the potential in the Schrodinger equa- 



tion ( [B.lp is positive for u < u = y 1 — 1/s 2 and negative for u > u . The solution 
to eq. ( |B.1| ), thus, decays exponentially in the interval (0, u ) and oscillates in the in- 
terval (uq, 1). Physically, the particle has to tunnel from u = to u = uq before it can 
reach the horizon u = 1. The imaginary part of G R is proportional to the tunneling 
probability, which is 

UQ 

Im G r (uj, k) ~ exp (^-2vo J du ^F(u)^j . (B.3) 

o 

For example, if to |fe|, then 

lmG R {uj,k) ~ e ~ a|fc|/T , w<|fe|, (B.4) 

where 

2r 

a = - v % ~ 0.835. (B.5) 

v^r (i) 

It is not immediately clear how to interpret this value of a. In a weakly coupled 
relativistic plasma, an external perturbation with spatial momentum (0, |fe|) can be 
absorbed by an particle with momentum (\k\/2, —\k\/2) so that its momentum becomes 
(| A; |/2, | A; |/2) and the particle remains on shell. The probability of such a process is 
suppressed by the Boltzmann factor of the particle, i.e., e~^^ 2T \ Since a > ~, there 
is some additional suppression in the strongly coupled CFT. 

The formulas presented above are confirmed by the more sophisticated Langer- 
Olver method, which is the WKB approximation for potentials with regular singular 
points (for more information and references see refs. j|, f45|1). 
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